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I. INTRODUCTION 

The Gregory-Laflamme(GL) instability is an s-wave instability for black p-brane solutions 
This instability is often quoted in the context of gravitational instabilities of black objects. For in- 
stance, the gravitational instability of ultra-rotating Myers-Perry(MP) black holes can be regarded 
as the one similar to the GL instability 3 1 . The nonlinear dynamical evolution of the GL instability 
has been numerically simulated in Ref. J] but as for the final state it is still under discussion. From 
view points of the membrane paradigm or field theory /gravity correspondence dual or 

analog fluid models of a black string and the GL instability have been attracting much attention 
in recent years [q- 

For the deeper understanding of gravitational instabilities and the duality with fluid models, 
it is useful to investigate responses to the variation or addition of physical parameters such as 
the electric charge, the number of dimensions and so on. One of the most fundamental and 
straightforward extension is to add rotation. In this paper, we consider the {D + l)-dim spacetime 
given by 

[D-dim MP black hole] x M , (1) 

and in particular focus on the I? = 4 case. We call these spacetimes MP black strings for D > A 
and the Kerr black strinq for D = A. 

Effects of rotation on the GL instability have been discussed by several authors 
The results of these previous studies suggest that, in general, rotation makes the GL instability 
stronger. We briefly summarize the known facts about the effects of rotation on the GL instability 
or instabilities on dual/analog fluid models. 



Ultra-spinning MP black strings 10|] 



When we consider the case of D > 5 and only one rotational axis, the angular momentum of 
the MP black hole is unbounded. For a very large angular momentum, the metric reduces 



to 



[D-dim MP blackhole] x M ~ [(£> - 2) -dim Schwraschild] x M'l (2) 



Since this is precisely the geometry considered by Gregory and Laflamme it is unstable. 
Therefore, we can expect that the GL instability persists even if rotation is added. 

Co-homogeneity-1 MP black strings 

MP black holes with all the angular momenta being equal in magnitude have the co- 
homogeneity-1 geometry (ie, the hypersurface of a constant radial coordinate is homo- 
geneous) in odd dimensions. Therefore, when we consider the geometry given by [co- 
homogeneity-1 MP]xM, it also has the co-homogeneity-1 symmetry in D dimensions as in 
the case of non-rotating black strings. We call these spacetimes co-homogeneity-1 MP black 



strings. In Ref. |l2l ]. from an analysis of stationary perturbations on co-homogeneity-1 MP 



black strings, it was predicted that the GL instability persists until extremality. In Ref. 0], 
the dispersion relation of the instability was numerically calculated and the prediction in 
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Ref. [Ij] was confirmed. The threshold wavelength and the instability time-scale become 
shorter as the rotation increases, that is, rotation makes the GL instability stronger for 
co-homogeneity- 1 MP black strings. 



Ray leigh- Plateau instability|14l| 



The Rayleigh-Plateau instability in a rotating fluid tube was discussed in Ref. [3] as a 
system dual to an MP black string. It is reported that the instability strength and threshold 
wavenumber increase as the rotation increases. 



Self gravitating cylinder [151 ] 



We can find a similar destabilizing effect of rotation on a self-gravitating cylinder in the 



4-dimensional Newtonian gravity |l5l ]. 



In this paper, we consider the Kerr black string with = 4. In this case, since Kerr black hole 
is not co-homogeneity-1, it is very difficult to solve full perturbation equations. Thus, we treat the 
rotation in a perturbative way by introducing a small parameter representing the rotation, that 
is, we consider only slow rotating black strings. We find that the rotation in the present case also 
makes the GL instability stronger. 

This paper is organized as follows. In SecHU we give a thermo-dynamical prediction comparing 
the entropy of a Kerr black string with that of a 5-dim MP black hole. Equations for the metric 
perturbations to second order in the rotation parameter are derived in Sec. IIIII and the results 
of the numerical calculation are presented in Sec. IIVI Some messy expressions are deferred to 
Appendices iB] and ICl A minimum review of spherical harmonics is given in Appendix |XJ We set 
the D + 1-dimensional gravitational constant and the speed of light equal to unity. 



II. THERMODYNAMICS 

First, we discuss the instability of stationary black objects from the thermo-dynamical point of 
view. We consider a Kerr black string with the length L along the z-direction and compare this with 
the 5-dimensional MP black hole which has the same total mass and the angular momentum as the 
Kerr black string. Since the horizon area is regarded as the entropy in black hole thermodynamics, 
if the horizon area of the Kerr black string is smaller than that of the MP black hole, the Kerr black 
string is thermodynamically less favorable than the MP black hole. Then, it is expected that the 
Kerr black string becomes unstable to a perturbation with the wavelength L along the z-direction. 

The metric of a Kerr black string is written in the form, 

,o A-a^sin^^,, „ . a ^r^ - A , , , 
ds^ = dr-2as\\i^e dtdcj) 

(r^ -ha2)2 - Aa^sin^e* o . , ,2 ^ , ^ ^ ,.2 , 2 
+- '—- sin^ edcl? + —dr^ + + dz^, (3) 

z-i ZA 

where < z < L with z = and z = L being identified, and 

S = r^ + a^cos^^, (4) 
A = - 2Mr + = {r - r+)(r - r_). (5) 
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The z =const., 4-dimensional section is given by the metric of a Kerr black hole. The tonal mass 
E and the total angular momentum J of the Kerr black string are given by 

E = LM, J = Ea. (6) 

The horizon area is given by 

where I and j are non-dimensional quantities defined by 

The metric of an MP black hole with single rotation can be written as 

= -dt^ + MML (df - aup sin^ edcP? + -^dr^ 
S Amp 

+ ^ (r^ + a^p) sin^ ^d<^^ + cos^ edu;^ (9) 

where 

Amp = r^ + a|ip - Mmp . (10) 

The total mass and the angular momentum are given by 

3 2 
E = -ttMmp , J = -Eaup ■ (11) 

The horizon area is given by 

tImp 2TT^MMpro 16 fs 9 



^3/2 - ^3/2 - 3 StT 4-^' ' ^^^^ 

where we have assumed Mmp > a|jp and vq is the horizon radius given by 

ro = Mmp - a^p . (13) 
We define the critical length Lcrit by the equality of the horizon area, that is, 

^BS = ^MP for L = Lcrit- (14) 
Then, fixing the total mass E, we find 

^BS > ^MP for L < Lcrit 

^BS < ^MP for L > Lcrit 

where 

Lcrit SVStT 



(15) 



_ V32 - 27pn =: krit- (16) 

VE Id 

Hence we expect that Kerr black strings are unstable against a perturbation with the wavelength 
L > Lcrit- The fact that the value of /crit is a decreasing function of j suggests that the critical 
wavelength of the GL instability becomes shorter, and hence the critical wavenumber //crit becomes 
larger if we add rotation. In the subsequent sections, by calculating the metric perturbations 
explicitly, we confirm this observation. 
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III. PERTURBATION EQUATIONS 

To avoid all possible complexities due to the coordinate singularity on the horizon, we use the 
metric in the Kerr-Schild coordinates given by 

. A-a^sm'^e^2 J 2asin2 0(r2 + a2- A) ^ ^ ^ ■2n^^ 
as = — dv + Mvar ^ avax — ^asm. Odxdr 

In Einstein gravity, the vacuum field equations are given by 

R^.u = ^■ (18) 

Here, we consider linear metric perturbations on the Kerr black string background. The metric 
may be decomposed as 

g = + A + • • • , (19) 

where 5(0) is the metric of the Kerr black string. 
Imposing the transverse traceless gauge condition, 

Vf,h^, = 0, h\ = Q, (20) 

we can rewrite the metric perturbation equations as 

Al/i^'^ := V^h^" + 2R^^f^k'^ = , (21) 

where Al is the so-called Lichnerowicz operator. 

As in the analysis of the original GL instability, we assume 

/j/^^ = , h^'' = Q. (22) 

Performing Fourier expansion along the symmetric spatial directions, we have 

f^ab ^ ^i{^,z+m4>)^^vj^ab^^^ (-23) 

where the subscripts a and h indicate the coordinates {v,r,6,x) and we have assumed that the 
frequency is pure imaginary as in Ref. 0]. Hereafter we concentrate on the m = case. The 
assumption of a pure imaginary frequency is correct for the non-rotating case [ll]. In our case, the 
validity of this assumption may be supported by the existence of such a solution. The equations 
of motion can be rewritten as 

{Ai^h)''^ = -fj,^h^^ + {Chy'' = ^ - i^h"^ + {ChT^ = (24) 

where 

~^ab ._ ^nv^ab (25) 
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and L is the Lichnerowicz operator in the Kerr spacetime. Gauge conditions are given by 

V,/i'^^ = 0, n = 0. (26) 

We solve these equations as an eigenvalue problem for /i with a fixed 0. 

As mentioned in Sec. HI we consider the slowly rotating case by introducing a small parameter 
e. We specify the deviation from the non-rotating case by 

M = M + e^M, (27) 
a = 2eM. (28) 

In the limit e — )• 0, the metric of a z =const. hyper-surface becomes the Schwarzschild metric. 
Since 

A = (r - 2M)(r - 2Afe^), (29) 

we have 

r+ = 2M, r_ = eV+. (30) 

In this way, the coordinate radius of the horizon is unchanged by adding the rotation. This fact 
makes our analysis much simpler because the coordinate value at the horizon boundary does not 
depend on the order of e. 

Fixing the value of we expand £, fi'^ and h as 

C = Co + eCi + e^Li + • • • , (31) 
= lil^ efij + eVi + • • • , (32) 
h = ho + ehi + e^h2 + ■■■ . (33) 

Since the value of fi should not depend on the direction of the rotation, we have fif = 0. Thus 
our goal is to get the value of /i2- The equations of motion can be schematically described at each 
order of e as 



e2 



-filhf + {Cohor'' = 0, (34) 
-/ig/if + (£0^1)'^' + {CihoT' = , (35) 
-filhf + (£0/^2)"' + {Cihir' + (£2/^0)"' - fJ^lK' = . (36) 



In the following subsections, we find that there exist a single master variable that satisfies a 
second-order ordinary differential equation at each order of e. The master variables must be care- 
fully chosen so that the equation can be easily solved by numerical integration. An inappropriate 
choice of master variables may cause 0/0 terms in numerical integration as reported in Ref. Q]- 
For this reason, we treat the = case separately from il. ^ cases. 
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A. = case 



Master Equations at order e" and 



As the same as the original work by Gregory and Laflamme, we consider only s-wave perturba- 
tions at order e". The metric perturbation can be written as 



tlQ 



(h^'"{r) /i^'"(r) 


y 










\ 






sin-^ d 



(37) 



Here we introduce iJ_ defined by 



(38) 



which is found to be a good master variable at order e*^. Using the gauge conditions (j26p and the 
equations of motion (j34p . we can derive the second-order ordinary differential equation for H- , 



Hi 



3/igr2 + 2 



+ 



4\ ^, ^ /igr {ulr^ + 4r - 8r+) 
r J ^ [r - r_|_) {fi^r^ + 2r — 2r+) 



(39) 



All the components of the metric perturbation can be expressed in terms of The explicit 

expressions are given in Appendix IB 1[ 

The metric perturbation at order is induced by the perturbation at order e*^ as given by 
Eq. (I35p . The only nontrivial induced components are found to be h^^ and h^^: 



( 



hf 



/i7(r)^ 
h\^(r) 


J 



(40) 



\/i7(r) h\''{r) 
Defining /(r) := h"^{r) as a master variable, we can derive the following equation: 



2^3 



2rj 



r'^{r — r+) 



-f + Si; 



(41) 



r+ {^rr\ {fi^r^ + 12) - Ar'^rl (^gr^ 7) + rV+ (/igr^ - 12) + Ur"^ - 12r^) , 

5*1 = z — ~ TTT^ — ?! — 7^ Z Z ^ H_ 



2r^{r — r^)^ (/^o'"^ + 2r — 2r^^ 
r+ {rl (10 - 3/x2r2) + 2rr\ (/xgr^ - ll) + r'^r^ (T/xgr^ + 6) - 2 {2iilr' + r^)) 

r'^(r — r+)'^ (/"o'^^ + 2r — 2r+) 



^. Master equation at order 

Equations at order are much more complicated than those to order e^. However, as we will 
see soon, we do not need to fully solve the perturbation equations to get the value of /U2. Instead 
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it is enough to focus on £ = modes in the spherical harmonics expansion. A minimal review of 
tensor and vector spherical harmonics is given in Appendix |3 In Eq. ()36p . the first and second 
terms correspond to the homogeneous part of the coupled differential equations. The operator of 
this part is the same as that at order e*'. This implies that different spherical harmonic modes 
of /i2 do not couple with each other within this homogeneous part. In addition, the final term 
proportional to has only the i = mode. Therefore, the value of /i^ can be obtained by solving 
only the £ = mode of Eq. (p6|) . 

Extracting the i = mode from Eq. (|36p . we obtain 



ab 



where 



"'2,0 



\ 




htoir) ■■= [hfir) 

which can be expressed as 

//i^^(r) hlf{r) 
hYir) hl'ir) 

\ 

We also extract the £ = mode from the gauge conditions. 

Here in analogy with the case of order e'', we introduce the following master variable: 



(42) 



(43) 



hl\r) 




(44) 



sin^ e I 



Ho 



1 



2hl 



(45) 



(46) 



9 



Then we obtain the master equation, 



r I — r r 



1 1 



/igr (/igr^ + 4r - 8r+) 
+ ^ ./no z z — r-n2- + o2 ; 



-f' + 



(r — r+) [fJ-^r^ + 2r — 2r^ 
4^Qr^(2r — r+) 



(47) 



+ 



;U^r3 + 2r - 2r+ (r - r+) (/Uq''^ + 2r - 2r+ 
1 / 



-/ 



3fiQr^{r — r+)(3r + r+) 



3r4(r - r+)2(/i2r-3 + 2r - 2r+)2 \ 
+/ior^(18r^ + 4r^r+ - 73rV^ + 14rVi]_ + 31rr+ - 12r^) 

+4r+(r - r+)(19r^ - 33rV+ - 3rVi + 17rri - 6ri)j + 6i4r^{2r - 3r+)(r - r+)^ ) iJ^ 



+ 



3/i2'^ (j^ — ^+) — 4rr+(/iQr + 6) 



+r^(/igr4 + 4^^r^ + 8) + 16r 



2r^3/i5r^ - 44) + rr^(-3Aior'' + Ufif^r^ + 208) 



-2rVJ|.(2/i^r^ + 35/igr2 - 24) + r^r+{-2fi^r^ + 23/i^r^ + 84^gr2 + 200) 
+3r5(/i[;r6 - S/z^r"^ - IG/igr^ - 16) - rV^ (/^e^e _ 26/i4r4 + 88^gr2 + 320) 



Exphcit expressions for the metric components in terms of H2~ are given in Appendix IB 31 



3. Asymptotic behavior 

In order to solve the master equations, we need to specify the boundary conditions. In the hmit 
r — )• oo, Eq. (p9|) can be approximated as 

Hi + ^H'_ - i^lH^ = . (48) 
Hence the general solution at r — )• oo is given by 

H- = Ci^— + C2—, (49) 

where we may assume > without loss of generality. Then the regularity requires C2 = 0. 
Thus decays exponentially at infinity, and so do all the components of the metric at order 
e'^. At order e, the homogeneous solution of the master equation is found to have the similar 
asymptotic behavior. Since the source term is given in terms of and H'_, this implies that the 
master variable f{r), and hence all the components of the metric perturbation at order e also decay 
exponentially at infinity. The same is true for the master equation at order since it has the same 
form as the master equation at order e". Thus all the components of the metric perturbations to 
order e"^ must decay exponentially at the infinity. 



10 



4- Near the horizon 

From the regularity at the horizon, we find H- = 0{{r — r^)^), f = 0{{r — r_|_)^) and H2- = 
0{{r — r^)^). Expanding the zeroth and first order master equations near the horizon, we obtain 
the following asymptotic solutions for i7_ and /: 



\3 



+0((r-r+)^), (50) 
f ^ ^ - ^+) + 12/3r+ + ^^y^ - 19/.gr^ + 94) ^ ^ ^ ^^^^ 

where a and /3 are constants to be determined. Without loss of generality, we may set a = 1 
because Eq. (j39|) is a homogeneous linear equation. The eigenvalue fiQ is determined so that i7_ 
decays exponentially at infinity, and the constant (3 is determined so that / decays exponentially 
at infinity. 

Since the zeroth order equation ()39p and the second order equation ()47p have the same form 
except for the presence of the source term ^2 for the latter, one can add the zeroth order solution 
to any particular solution H2- of Eq. ()T7|) . Using this freedom, we can set H2-{r+) = by 
adding a term proportional to appropriately. Then expanding Eq. (j47p near the horizon, we 
find 

H2- = (r - r+) (^4/3 + r+(^^+ - 1^J + Q {{r - r+f) . (52) 
As before, the value of ii\ is determined so that H2- exponentially decays at infinity. 



B. case 

In the O 7^ case if we would use the same master variables as those in the $7 = case, we 
would suffer from 0/0 terms in the master equations at order e'^ and e^, as pointed out in Ref. [2I. 
To avoid numerical errors due to these 0/0 terms, as in Ref. Q], we rewrite each master equation 
at order e*^ and as two coupled first-order differential equations by introducing an additional 
variable. 

At order e'', we define a function H as 

H := -hl^ + -X^K- (53) 
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Then, we obtain the following equations for and H: 



r 



n {filr^ + r+) 

(^gr3(3r+ - 2r(J]r + 1)) + 4f)2r4 - 2rr+(J7r + 3) + 6r2 ) 



r 



2 



4(r — r+)2 (/^o'^^ + 1"+) 
{fj,lr^{2r{nr + 1) + r+) + AQ'^r^ + 2rr+(ilr + 2) - 2r^) ^ ^^^^ 
2r(r — r+) [^qT^ + r+) 

These equations are free from 0/0 terms. The components of the metric perturbation are expressed 
in terms of H- and H as shown in Appendix O 

At order e^, we can use the same master variable as the one in the i7 = case, / := h^^. It 
satisfies 

_^// ^ ( ^ 2(l]r + 1) _ 3\ ^, ^ fiy-2filrr+-24n^r^-12nr+ ^ 

V/iQr + 4il r+ — r r/ r2(r — r+) (^gr + 4^7) 

+AH{r; r+, ^0, ^^)^^ + {r; r+, ^o, !^)^^- , (56) 

where explicit expressions of Anir; r^, fiQ, Q) and Ah_ {r; r+, /j^q, 17) are given in AppendixO The 
other component of the metric perturbation is also expressed in terms of / in Appendix ICl 
At order e^, we define H2 as 

1 

1 " 

Then, we can derive 



, {An'^r^ + /ig(3r+ - 2r(0r + l))r^ - 2(J7r + 3)r+r + 6r^) 

i^2— ~ 7r~/ w — \ ^2- 

2r(r — r+j(//Qr'^ + r+j 



H „ , 9 o ^ -n2 



r(r//Q + 4f])(/iQr3 + r+) r'^{r^'^ + 4il)(/iQr3 + r_|_) 

+BH{r; r+, /xo, ^2, f^)-f^ + (?^; ?^+, fJ.2,^)H-, (58) 

, r2r(4J7^r'' + 4//o(?' - r^)r^ - r\) 

H2 = —} 757 — 9-^5 ^ -"2- 

4(r — r+j^(/UQr'^ + r+j 
(4J7V^ + /ig(2r(17r + 1) + r+)r3 + 2{VLr + 2)r+r - 2r^) ^ 
2r(r — r+)(/XQr'^ + r+) ^ 
_^ 2VL{2{riil + J])r2 + r+)r2 ^ 2Sl(2(r/xg + 2f](r (r^g + Vt) + 2))r^ - 40r+r2 + r\)rl ^ 

r{rjiQ + 4r2)(/XQr3 + r+) r^{rfj,l + 4r2)(r — r+)(/iQr3 + r_|_) 

+C/^(r; r+, /xo, /U2, 0)i? + C/f_ (r; r+, /xq, /U2, ri)i?_ , (59) 

where explicit expressions of BH{r;r-^, hq, fj,2,^), BH_{r;r+, fiQ, jj,2,^), Ci^(r; r+, /xq, Ai2, f^) and 
Cffir', r+, Ho, /^2, ^) are given in Appendix [Cl The components of the metric perturbation are also 
shown in AppendixO 
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At infinity, all these variables must exponentially decay. We can derive asymptotic solutions 
near the horizon as in the = case. We describe these asymptotic solutions near the horizon in 
Appendix [Cl 

IV. RESULTS AND DISCUSSION 

Since a second-order ordinary differential equation can be regarded as a set of two coupled 
first-order differential equations, we have six coupled first-order ordinary differential equations for 
each value of We have to determine the values oi ^q, (3 and ^2 so that the metric perturbation 
is regular at horizon and decays exponentially at infinity. It can be done numerically by a shooting 
method. We have solved the ordinary differential equations from the horizon to infinity using the 
4-th order Runge-Kutta method and searched for appropriate values of /xq, /3 and ^2 for each value 

of n. 

Here, we show the numerical results by fixing the mass parameter M of the Kerr black hole 
at the 4-dimensional section, although the numerical integration is done by fixing the value of r_|_. 
That is, we plot Q,M as a function of fJ,M, where /i = \/fJ'Q + /^i^^' ™ ^^S- E where 

QM = nM{l + e'^), (60) 
= ^^2 + ^2,2^(1 + e') = (1 + (1 + e') + 0{e') . (61) 

It is clear from the figure that both the critical wavelength and the growth time-scale of the 
instability become shorter as the rotation is increased. In other words, rotation makes the black 
string system more unstable. This result is consistent with previous works in similar situations as 
listed in Sec. HI 

Although our results are consistent with all the previous results, an intuitive understanding 
of our results seems rather unclear. If we fix the mass parameter of the Kerr black hole in the 
4-dimensional section, the horizon area decreases as the angular momentum is increased. In this 
sense, the typical time-scale or length-scale becomes shorter. This is consistent with our and 
previous results, and this was in fact what we obtained from thermodynamic arguments. However, 
in contrast, if we focus on the surface gravity, it decreases as the angular momentum is increased. 
In this respect, the typical time-scale or length-scale may become longer. Actually, when we add a 
charge to a neutral black brane solution, it is known that the GL instability gets weaker Qj- This 
is totally opposite to the behavior we obtained. Apparently we need a better physical picture of 
the GL instability with rotation. 
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Appendix A: Spherical Harmonics 



Full expressions for the tensor spherical harmonics can be found in Refs. 16Hl8t|. Here, we give 
a minimum review necessary for the present work. We represent tensors by using the components 
of the coordinate system (v,r,9,(j)), where v is an advanced null coordinate. First, we introduce 
the vector spherical harmonics. We decompose any covariant vector field v as 



V = ^ [Aimaim + Bimblm + ClmClm + Dimdlm] , 



(Al) 



Lm 
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where A, B, C, D are the coefficients for bases a, h, c, d given by 

aim = -i(l]m,->lm,0,0), (A2) 
him = (0,yi„,0,0), (A3) 

Cim = -j=^={Q,Q,dBYim,d^Yim) , (A4) 

1 1 

dim = —j===={Q,Q,^^d^Yim,-sm.edeYim). (A5) 
These bases are orthonormal in the inner product 

{v,u) := / t^aUbdn, (A6) 



where r] is the Minkowski metric. 

We focus on the I = m = mode. In this case, we have the following two non-zero bases: 

a := aoo = ^{l, - 1, 0, 0) , b := boo = 77^(0, 1, 0, 0). (A7) 
Thus, we can extract the I = m = mode from a covariant vector v as 

[v]i=o = {a,v)a+{b,v)b. (A8) 

In the same way as in the vector spherical harmonics, we can decompose a rank-2 covariant 
tensor v as follows: 

V = ^ I^^Zma/m + Bimhlm + Clm^lm + Dimdim + Eim^lm + Flm^lm 
Lm 



~^GimSlm ~l" Himii-lm ~l" IlmXlva ^ Jlnii 



lm 



(A9) 



where a, b, c, d, e, f, g, h, i, j are bases of the tensor spherical harmonics and 
j4, C, • • • , /, J are the coefficients for these bases. Explicit forms of the bases are 
given in Ref. [18]. These bases are orthogonal in the inner product 



< V, u >:= / r^'^'^rlx^VabUcddn. (AlO) 
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For m = Z = 0, we have the following four non-zero bases: 

/ 1 -1 o\ 



a := ao ,0 



b := b, 



'0,0 



c := co,o 



do,o 



1 



20F 



2V27r 



20F 



2V27r 



-110 


\ / 

/ -1 o\ 

-12 

y 

/ \ 

10 



\ J 

/ ^ 


1 
\^ sin^ e J 



We can extract the I = m = mode from v as 

[v]i=o =< a, V > a+ < b, V > b+ < c, V > c+ < d, V >d. 



(All) 



(A12) 



(A13) 



(A14) 



(A15) 



Appendix B: metric perturbation for f2 = 

Here we give explicit expressions for the components of the metric perturbation at each order 
of e in the case = 0. For notational simplicity, e is set to unity in the following expressions. 

1. Metric perturbation at order e'^ 



riQ 



hi' 



r^(2r — 3r_( 



-H' 



{fi^r^ + 4r - 8r+) 



(r — r+) {hqT^ + 2r — 2r+) (r — r+) {fiQr^ + 2i — 2r+) 



i7_ 



r{2r — 3r+) 
2;u§r^ + 4r - 4r- 



Sr. 



fiff^ + 2r - 2r4 



-H. 



(2r-3r+)(r-r+)^, ^ (r - 3r+)(r - r+) ^ 



2/XQr3 + 4r - 4r+ 
1 



r (|ugr3 + 2r - 2r+) 



2r2 



(Bl) 
(B2) 
(B3) 
(B4) 
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2. Metric perturbation at order 

h7 = H' + — M:^^!±) (B5) 

^ 2/igr4 + 4r2 - 4rr+ (//^r^ + 2r - 2r+) 

3. Metric perturbation at order 

^rr ^ (2r-3r+)(r-r+) ^, ^ (r - 3r+)(r - r+) ^ ^ 
^ 2/ugr3 + 4r - 4r+ ^" r (/7,^r3 + 2r - 2r+) 

^ 2r^(r-r+)2 ^ 2r^(r - r+)(2r - r+) 

r (/igr3 + 2r - 2r+) ^ (/xg^s + 2r - 2r+) ^ 

+ |^r+ (-2lM^r^ + 31fiyr+ + 2r'{^ilrl - 3) + 

r^+iAO - 9//gr^) - r2r^(3//^r^ + 59) + 4rVi[ + 9rr^ + 6r^) 

-3;z^r^(2r - 3r+)(r - r+) j 6.3(^2,3 + 2r - 2r+)^ 

+ ^r+(^-36/x^r^ + 79nlr^r+ - 4r^(4/xgr2 + 9) 

+8rV+(14 - Sulrl) + rV^(3/xgr^ - 176) + 34rV^ + 48rr| - 6r%J 
-'"'^'•°"--^'-^'<'----^>) 6r^(^g..f2r-2r,F - C^"' 



r2(2r-3r+) _^ r (^ugr^ + 4r - 8r+) 

(r — r+) (/Ugr^ + 2r — 2r+) ^ (r — r+) (/Xgr^ + 2r — 2r+) 



4r r 2 _ ^ 4r^(2r-r+) 



IJ,Qr^ + 2r - 2r+ (r - r+) (//q?'^ + 2r - 2r+) 

+ |^r+(-15/xgr^-3r4(/x§r2 + 3) 

+2/r+(ll//^r2 + 5) - rV^(4/x^r2 + 53) + 6r^ + 34rV^ + 6r^) 
— 3//|r^(2r — 3r+)(r — 



3r(r - r+)2(;ugr3 + 2r - 2r+)2 
+ 1 r+ ( 3rr|(Aigr2 + 4) + 2r2ri|.(41 - 3nlr^) 



-6r^(5Aigr2 + 4) + rV+(6lMor^ + 64) - 2rV2 (ll/x^^^ ^ _ g^5 ^ 

-6;.ir^(r - 3r+)(r - ^+)) 3,2(, _ ,^)2(^?;3 + 2r - 2r+)2 ' ^^^^ 



r(2r - 3r+) ^, ^ (r - 3r+) ^ ^ 
2//gr3 + 4r - 4r+ ^" /Ligr^ + 2r - 2r+ ^~ 
^ 2r^(r-r+) ^ 2r^(2r-r+) 

^2^3 + 2r - 2r+ ^ /n^r^ + 2r'^ - 2rr+ ^ 

+ - 3rr|(/i^r2 + 3) + 2rV+(ll/i^r2 + 5) - rV^ (4/xgr2 + 53) + 6r^ 

+ |^r+ (3rr^(/igr2 + 4) + 2rV^(41 - 3/x§r2) 

-6r^(5/x^r2 + 4) + r\+{61nlr'^ + 64) - 2rV^(ll/x^r2 + 76) - 6r^) 
— 6/x|r^(r — 3r+)(r — r+) 



6r^(r — r+)(/XQr^ + 2r — 2r+)^ ' 



h^O ^ _ ^ ^/^«x _ r-+(2r - 3r+)(3r2 - 2rr+ + r^) ^, 

^ 2r2 3r3 ^ 6r^{r — r^){fiQr^ + 2r — 2r^) ~ 

+r+ {-3fj,y - 2fiyr+ + 2r^ {filrl - 6) + r\+{filrl + 24) 

-6rVi + 4rri - 2ri') — , ^ . 

+ + +/6r5(r-r+) (/x^r3 + 2r-2r+) 



Appendix C: expressions for Q ^ 
1. Metric perturbation at order e° 



r [2ij.lr^{r — r+) + 4ri^r^ + r+(r+ — 2r)) ^ 
2(r - r+)2 (|U^r3 + r+) 

r (//2r(2/-(Q/- - 1) + 3;-+) + 29. {29.r^ + r+)) 



n{r - r+) {ulr^ + r+) 



_ 40V + r+(3r+ -4r) ^ _ /xgr(2r(nr - 1) + 3r+) + 2Q(2^^r^ + r+) ^ 

° ~ 4(r - r+)(/igr3 + r+) ~ 2ri(^gr3 + r+) ' 

^„ _ (r - r+) (/ig(2r - 3r+) - AQ'^r) ^ (40V^ - 4rr+ + 3r j) ^ 



1 
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2. Metric perturbation at order e 



+2//gQ(^2J^r'^ + r+ (-4r2 + 8rr+ - 3r^)) + 4f7V+(2r2r + l)(2r - r+) j 

+2nr+(^-2r^(fir(4nr + 1) - l) + 2rV (^^(2^^^ + 1) + 3) - Urrl + Qr\^ H _ 

+ /x2r + 4S^^i r(/x2r + 40) " ^^^^ 
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3. Metric perturbation at order 



h7 



2(r - r+)(2r{rnl + n{i4r'^ + 2nr + 4)) - (r/x^ + 60)r+)i 
r(r/xg + 4J^)(/Ltgr3 + r+) ^ 



402r^ - 4r+r + 3r% 

o — A ^2 

Afi^r^ + 4r+r 



(r - r+) ^(3r+ - 2r)/x§ + 4J^Vj 
i^s>- . _ / „ „ ^ H2 



2n {fiy + r+) 



6iJ.l^l{riJ,l + 40) (r - r+)(4J^V + r+(3r+ - 4r))j 



+r+ r^(r - r+) (^24J^ (J^r(J7r - 1) + l)r^ + 4(Or(J^r(2J^r + 3) - 15) + 6)r+r^ 

-2(Or(6Qr - 17) + 15) r^r^ - 3(20r + 9)r^r + 3(9 - 160r)4)/x^ 

+J]r (^24J](Jlr(2Jlr(J]r + 3) - 3) + 4)r'^ + 16(J]r(Jlr(2Jlr(J]r - 2) + 3) - 24) + 6)r+r^ 
+2(2Jlr(4Jlr(J7r(Jlr - 5) - 8) + 79) - 105)r^r^ - 4(Jlr(1917r + 16) - 42)rlr^ 
+3(4Qr(7Jlr - 1) + 39)r|r2 + 3(16J^r - 81)r^r + 102r^)jU^ + 4^^^(24(72 (2J^r(J^r + 1) + 1)^ 
+16n{2n^r^ - 6nr - 3)r+r^ + 2(2J7r(4nr(Qr(J7r - 3) - 2) + 9) + 3)rir^ 



+4(17r(6 - 7J)r) + 9)riv^ + 3(2f]r - l)(6f]r + l)rV^ - 51r3.r + 30r 



12ri2r3(r/xg + 4J7)(/x2r3 + r+f 



6/x^r2(rAt5 + 4Q)(r - r+){AQ^r* + (2r - 3r+)r+) 



3r^(r - r+)(lOnr^ - 2{D.r+ + 2)r^ + 2r+(4 - 3nr+)r'^ + r^(3 - 2fir+)r - 9ri|.)//^ 

+0r (l2(2n3r^ - 14fir + 3)r^ + 4(fir(Or(6J]r + 5) + 48) - 30)r+r'' 
+2(J]r(2J]r + 21) + 50)r^r^ + (83 - 102nr)rlr^ - 3(8Jlr + 59)r|r + 48r^)/x^ 
+2J^2 (^12(fir(2J^r(3fJr + 1) - 7) - 2)r^ + 4(fir(J^r(16fJr + 13) + 18) + 12)r+r^ 
+2(Or(2Qr(20r - 1) + 13) - 14) r^r^ - (176nr + 5)r|r2 + 3(10Qr - 27)r|r + 24r^)/x§ 

+8n^(l2n{2nr{nr + l) + + 4f)(r2r(41]r + 3) - 6)r+r^ 



+2(Or(2Qr(20r - 3) - 5) - 6)r^r2 - (20nr + 3)r^r + 3(2J^r - l)r^ 



(C6) 
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(r/x2 + 4J7)(//2r-3 + r-+) ^ (r/xg + 4J^)(r - r+)(/x2r3 + r+) ^ 



r(4ri2r* + 2/XQ(r - r+)v^ + r+(r+ — 2r)) ^ 
2{r — r^yip^r^ + r^) ^ 
I ^ |^ r(^g(2r-3r+)-4^^V) _ 3 | 



r — r+ 



0(//gr3 + r+) 
1 



X [ 6filn{rfil + 4n)(r - r+)(4QV + r+(3r+ - 4r))r^ 



+r+ r^{r - r+) (2m{nr{nr - 1) + l)r^ + 4(Or(J^r(2J^r + 3) - 15) + Q)r+r^ 

-2{nr{6nr - 17) + 15)r^r2 - 3(2J^r + 9)r^r + 3(9 - 16J^r)r^)//^ 

+$7r(24f7(r2r(2r2r(Jlr + 3) - 3) + 4)r'^ + 16(Jlr(J]r(2r2r(f7r - 2) + 3) - 24) + 6)r+r^ 
+2{2nr{Anr{nr{nr - 5) - 8) + 79) - 105)r^r'^ - 4(Jlr(19J]r + 16) - 42)ri|.r^ 
+3(4Jlr(717r - 1) + 39)r^r^ + 3(16Jlr - 81)r^r + 102r'j:)^g + 4n'^{2An^{2nr{nr + 1) + l)r^ 
+16fi(20V^ - enr - 3)r+r^ + 2(2J^r(4J^r(J7r(J^r - 3) - 2) + 9) + 3)r^r'' 

+4(fJr(6 - 70r) + 9)r^r^ + 3(2fir - l)(6^^r + l)4r2 - 51r^r + 30r^) 1 iJ_ 



6n'^r{rnl + 40) (r - r+)2(//gr3 + r+)^ 
X ( 6/xiO(r/x^ + 40)(r - r+)(4J)V'^ + (2r - 3r+)r+)r3 



+r+(-3r^(r - r+)(10^^r^ - 2(nr+ + 2)r^ + 2r+(4 - 3^^r+)r^ + r^(3 - 2nr+)r - 9r^)/x^ 
+nr(12(20V^ - 14J^r + 3)r^ + 4(J^r(Or(60r + 5) + 48) - 30)r+r^ 
+2{nr{2nr + 21) + 50)r^r^ + (83 - 102Or)r|r2 - 3(8fJr + 59)4r + 48r%)n^ 
+2n^{12{nr{2nr{3nr + 1) - 7) - 2)r^ + 4(Or(J^r(16J^r + 13) + 18) + 12)r+r^ 
+2{nr{2nr{2nr - 1) + 13) - 14)r^r3 - (176J^r + 5)rlr'^ + 3(10Or - 27)r|r + 24r^)//^ 
+8J7^(12J7(2f]r(Jlr + 1) + l)r^ + 4J](Jlr(4f]r + 3) - 6)r+r^ 



+2(fir(2fir(2fir - 3) - 5) - 6)rlr^ - {20nr + 3)rlr + 3(2J^r - l)r|)) H, 



(C7) 
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2(r^2 ^ 2n){r - r+)rl , 4r(r/xg + ^{iilr^ + 2^r + 4))r2 - 2(rMg + 6^^)r3 



(r^2 + 4J7)(^2^3 + 



1 / r(^g(2r-3r+ 
2V 0(M^r3 + r+) 



r(r^o + 4r2)(//5r3 + r+) 

40V + r+(3r+-4r) ^ 
Z I tl2 H 77 , , 9 „ r-"2- 

4(r - r+)(//5r'^ + r+) 



X [ 6/xiO(r/x^ + 40)(r - r+)(4J)V^ + r+(3r+ - 4.r))r^ 



+r+ r\r - r+) (^240(J)r(Or - 1) + l)r^ + 4(Or(J^r(2J^r + 3) - 15) + 6)r+r^ 

-2{nr{6nr - 17) + 15)r^r2 - 3(2J^r + 9)rlr + 3(9 - 16J^r)r^);U^ 
+J]r(^24Jl(J]r(2J7r(Jlr + 3) - 3) + 4)r'^ + 16(J^r(Jlr(2Jlr(J]r - 2) + 3) - 24) + 6)r+r^ 
+2(2J^r(4J^r(fJr(J^r - 5) - 8) + 79) - 105)r+/ - 4(J)r(19fir + 16) - 42)r^r^ 
+3(40r(7Qr - 1) + 39)r^r2 + 3(160r - 81)r|r + 102r^)/ig 

+40^(^2402 (20r(J^r + 1) + l)r^ + 16J^(20V^ - 6J^r - 3)r+r^ 
+2(2f]r(4fir(Or(f]r - 3) - 2) + 9) + 3)r^r^ + 4(J^r(6 - 7J]r) + 9)rij_r^ 

+3(2^^r - l)(60r + l)4r2 - 51r^r + 30r^) ^ 



X ( 6i^ln(ri^l + 4f])(r - r+)(40V^ + (2r - 3r+)r+)r^ 



-3r^{r - r+) (^lOJ^/ - 2(0r+ + 2)r^ 

+2r+(4 - 3nr+)r2 + r^(3 - 2J^r+)r - 9ri|.) fi^ 

+nr(l2{2n^r^ - UQr + 3)r^ + 4(fir(Qr(60r + 5) + 48) - 30)r+r^ 

+2{nr{2nr + 21) + 50)rlr^ + (83 - 102nr)rlr^ - S{8nr + 59)4r + 48r^) /x^ 

+2n'^(l2{nr{2nr{3nr + 1) - 7) - 2)r^ + 4(J]r(Qr(160r + 13) + 18) + 12)r+r'' 

+2{nr{2nr{2nr - 1) + 13) - 14)r^r^ 

-{nenr + 5)rlr'^ + 3(10Jlr - 27)r^r + 24r^) //g 

+8^^^ (|l2J^(2J^r(fir + 1) + l)r^ + An{nr{4nr + 3) - 6)r+r^ 
+2(r2r(2fir(2fir - 3) - 5) - 6)rlr'^ - (20Or + 3)rij_r + 3(2fir - 1)^) j . 



(C8) 
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+ ^± 

6J]r4(r - r+){nlr + 4n){iilr^ + r+) 

X ^/Ltor^ (6r^{nr - 1) - 2rr+(i7r + 4) + 13rV+ + 3r^) 

+2//gJ7 (^6r^(3J7r - 2) + rV+(4f]r + 21) - 2r'^rl{3nr + 2) + 9rrij_ - 6r^^) 

+8^^2(6^^r^ - 2r|(^^r + 1) + 2rr^(nr + 1) + 3rV+) j fl" 

+ ^± 

12r5(r - r+)2(/xgr + 4:n){nlr^ + r+) 

X ^—2iJ,Qr^{r — r+)(r + r+)(3r^ — rr+ — r^) 

+^^r(-12J^r^(J^r + 2) + 2r-^r+(8J7r + 3) + rV^(4fir(rir + 6) - 15) + 2rV^(9 - lOJ^r) 



+rrX{4nr - 5) - 2r^) + 4Q(^-120V - rV^(2i7r(20r + 1) + 5) 

+2rV^(J^r(2fJr + 1) + 5) + 6rV+ - 12rr| + 4r^)^ (C9) 



4. Ah, Sh, Bh_, Ch and C^. 



Aij(r;r+,/xo,^^) 

2J]r3(r - r+)3 (^g/^ 4^^) (^2^3 + ^^) (/^o^'^+(^ " " ^+) 

+2/x^r(-2r^(J^r(J7r + 1) - 3) - r^r+{nr{4nr + 1) + 6) 
+2rV^(J^r - 7){Vtr + 1) + 9rri|.(0r + 2) - 6r|) 
-8^^0(r^(0r + l)(3J^r - 2) + r3r+(J]r(6J^r + 11) - 3) 
+rV^(19 - Qr(3J7r + 2)) - rr\{9.r + 14) + 3r^) 

{-2^r^{^r + 1) - 4J^r\+(Qr + 2) + rjj.(3 - 2J^r) + rr^(2J^r(J^r + 3) - 5))) , (CIO) 
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AH_{r;r+,iio,Q) := 

+2rV+(317r+ + 13) + 17rr^ - 21rJ|_) 

+Afilr{2nr'^{nr - l){nr + 6) + r'^r+(Jlr(^7r(4J7r + 1) + 34) - 10) 
~2r^rl{nr{nr{nr + 3) + 8) - 6) - rV!|_(f^r(J7r + 19) - 10) + 2rr|(6fir - 11) + 9r^) 
+8n{4n'^r'^{nr + 1) + 2r\+{nr{nr{'inr + 7) - 1) - 5) 
-2r^rl{nr + l)(2J]r(Jlr + 3) - 3) + 3rVi^(20r(f^r + 2) + 7) 

-r4(6fir + 31) + 124)), (Cll) 

-Bi?(r;r+,/xo,/X2,0) := 

o rr-T-^ ^i^lA^irfil + 4^^)(r - r+)(/x^r6 + 4J^V^ 

+2nlr+r^ + 2{r - r+)r+)r^ + r+(r*^(r - r+)(6(rir + l)r^ + (2Jlr - 3)r+r - 9r+)/io 
+r''^(6(f7r(2f7r($7r + 2) - 1) + 2)r'' + 2{nr + 3)(2f7r(2Jlr - 1) - 5)r+r^ 
+ (rh-(4f7r(r2r + 1) + 21) + 6)rlr^ + (43Jlr + 30)rJ^r - 18(r2r + l)r:|)//o 
+4f]r(3(8Jl^r-^ - UUr + 3)r^ + 4(Jlr(J]r(5J]r + 16) + 8) - 6)r+r'^ 

+ {2nr{nr{2nr + 11) + 30) + 9)rlr^ - 2(Jlr - 2)(f]r + 6)ri|:r2 - 6(317r + 5)r|r + 9r^)//^ 
+4Jl2r(6(Jlr(2Jlr(5Jlr + 1) - 7) - 2)r'^ + 4(f^r(2J7r + 3)(7J]r + 5) + 6)r+r^ 
+(Jlr(217r + 9)(2Jlr + 11) - 36)rlr^ + (42 - nr{8i^r + 87))ri|_r + 12{i^r - 2)r|)At^ 

+i6n^{6n{2nr{nr + l) + l)r^ 

+4Q(fir(2fir + 3) - 3)r+/ + 0(40^^ + 15)r^r2 - 3(0r - 2)r^r - 6rj))), (C12) 
BH_{r;r+,HQ,ii2,^) ■= 

l2nr\r^.l + 40)(r - r+)^(^gr3 + r^)^ (^+(-3/^o(^ - -+)(8^^^(- + 
+2J7(-3r^ + 3r+r^ + 3r+r + r\)r + (4r - 3r+)r+(r + r+))r^ 

+/i^(-24J](J7r(J7r(J^r + 5) - 4) + l)r^ - 4(J^r(J]r(J7r(6J^r + 13) + 54) - 18) + 6)r+r^ 
+2(17r(3Jlr(14J^r + 5) - 31) + 21)r+r^ + {VLr{2VLr{%Q.r + 39) + 57) + Q)r\r'^ 
+(^^r(240r - 1) - 42)r|r + 6(3 - 4fir)r^)r^ + 2iilQ.{-l2Q.{?,ilr{2ilr{9.r + 1) - 1) + 4)r^ 
-A{ilr{Q.r{^r{lQQ.r + 29) + 24) - 66) + 12)r+r^ 

-2{nr{Q.r{2Q.r{2Q.r + 3) - 21) + 86) - 45)r^r'^ + {9.r{lQQ.r{2Q.r + 7) - 27) - b{))r\r^ 
+{Q.r{l - MQ.r) - 50)r^r^ + 2(12Jlr + 47)r^r - 42r^)r + m'^{-l2if{2ilr{Q.r + 1) + l)r^ 
+AQ.{nr{l2 - nr{AQ.r + 3)) + 6)r+r^ - 2{Qr{Q.r{2Vtr{2nr - 3) + 9) + 12) - 3)r+/ 
+{Q.r{AQ.r - 9) - 2)r%r^ + (5r2r(3 - 2Jlr) - 14)r|r^ + 16r^r - 6r^)) 

-Q^JLl^r^{r^ll + 4Q)(r - r+)(4QV^ + r+(3r+ - 4r))), (C13) 
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CH{r;r+,iJ.o,iJ.2,Cl) := 

+r+(-r^(r - r+){<d{7Vtr - 4)r^ - 2{9.r - 15)r+r^ + (51 - 10J]r)r^r - 3(2Jlr + 21)ri|_)//^ 
+r2(12^](2^]'V^ - 20Jlr + 7)r^ + 4(Jlr(J7r(f]r(6Jlr + 5) + 56) - 54) + 3)r+/ 
+2(J]r(Jlr(2Jlr + 1) + 34) - '^)r\r^ + (J]r(243 - 70J7r) - 36)rJ|_r2 

+ (78 - 9.r{2AQ.r + 353))r^^r + 12(1017r - ?>)r%)i4i + 20(12fi(Or(20r(30r + 1) - 11) - 2)r^ 
+4(^]r(Jlr(J]r(16Jlr + 13) + 18) - 3) + 6)r+r^ 

+2{Q.r{9.r{29.r{29.r - 1) + 37) + 2) - 21)r^r'' - {Q.r{272Vtr + 81) + 2Q)r%r'^ 

+ (Or(78J]r - 25) + 82)r^^r^ + 4(617r - 23)r!|_r + ?>Qr\)iil + 8J72(12J]^(2J7r(Jlr + 1) + l)r^ 

+4Q(rir(Jlr(4Jlr + 3) - 6) - 3)r+r^ + 2{Q.r{ilr{29.r{2nr - 3) - 5) - 6) - 3)r^r^ 

+(J^r(9 - 20fir) + 4)r|r2 + (3fir(2i^r - 9) - 14)4r + 2{mr + h)r%))), (C14) 

Ci?_(?';?'+,M,/^2,f^) := 

o o o r:r(r+(-12/x^(2r - 3r+)(r - r+)2(r + r+)r^ 

24r4(r/xg + 4r^)(r-r+)3(/i2r3_|_^_^)2V +v /^uv -h;v +j \ +j 

-nl{r - r+){24n{^r{nr - 2) + 5)r^ + 4(Or(Or(20r - 3) - 51) + 24)r+r^ 
-~2{VLr{<oVLr + 7) + 48)r^r2 + 3(460r - 47)r|r + 15(9 - mr)rX)r^ 

+/x^(-2402(Or + 4)(2fir(J^r - 1) + l)r^ - 8^^(^^r(2^^r(2^^r(Or - 2) + 3) - 51) + 48)r+r^ 
+2(Or(2^r(40r(^^r(5 - Or) + 2) - 79) + 381) - 24)r^r'' 

+2(20r(J^r(670r + 22) - 66) + 39)r|r^ - 3(3J^r(40r(50r + 1) + 85) - 20)r|r2 
-3(J^r(16Qr - 305) + 54)r^r + 6(12 - hZQ.r)r%)r'^ + 4^^(-24^^3(2^^r(^r + 1) + l)r^° 
+8^^2(-4^^V^ + 12Qr + 9)r+r^ - 2(^^r(2^^r(4^^r(^^r(^^r - 3) - 2) + 9) + 3) + 24)r^r^ 
+30(40r(5 - 30r) + 5)r|r'' + 2(20r(J^r(70r - 12) - 9) + 39)r|r^ 
-3(J^r(8nr - 13) + 38)r^r2 + 6(19 - h^r)r\r - 36r^)) 

-6/i^0r^(r/x§ + 4^^)(r - r+)(4QV'' + r+(3r+ - 4r))). (C15) 
5. Expansion of meister variables near horizon for f2 7^ 



Q.r+{2VLr+ - 1) ^ nl^rX + VL 



(r - r+) (2^2^+ + 4^?) 2^2^+ + 4J1 
fi(r - r+) (-2^7r+ (^^^j - 7) - ?,^lr\ + 8^gr2 + 402^2 + 4) 



= 1 + 



8r+(J]r+ + l)(2f^r+ + 1) (^2^+ + 2^) 
(r - r+) {^ilr% - 4^gr+ - 8J^2r+ - QQ) 
r+{2rtr+ + 1) {nlr+ + 2^) 



+ 0((r-r+)2), (C16) 



+lilr+ {tilrX - ^ixlrl + 32) + 481^^2 ) + 0{{r - r+f), (C17) 
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1 ^ ^ 4r2 (f]r+ + l)2(217r+ + l)(2r?r+ + 3)(/i2r+ + 2f))(^2^+ + A^l) 

{Sp^Mln^rl - 4J]r+(4J7(9/3 + + 109/Uo) - 2r^(2/UoJl(24/3 + 209O) 

+803(51^5 + + 47/i^) + 8/3^^J^V^(3/ug + 21)2) _ (/3(6^^ + 248^gJ]2 ^ gigj^^) 

+53/i^O + 210^gj]^ + ll2if) + ^llr%{2^lla{ll|5 + b^) + 8/igj^^(4/3 + 0) - 192/3J]5 - 

+r^(-2/i^J7(16/3 + 317) - 24^gf]3(3g^ ^ gj^) ^ - 21/3) + 19/i^) 

+r^(/3(6/i[j - 4/i^f]2 _ 704^2 _ 3345^6) ^ ^^g^fSj^ ^ 24^^f]3 + 32^gJ7^) - 408Jl2) 

+0((r-r+)2), (C18) 



(r — r+) 



3r+(f]r+ + l)(2J7r+ + 1)2(^2^+ + 2J7)2(/i2r+ + 4f]) 

(Aigr3 (12/3r+ + ^^^2 + 3^2^2 _ 21) _ SJ^V^ (-24/3r+(3/i2r2 + 4) + 2^ilrX - 21^2^2 

-24//|r2 - 62) - 4fiV+(-24/3r+(3^^r^ + 12^gr2 + 5) + + 6^^r^ 

-/^gr2 + 125) - eO/^^r^ - 64) + 2Q^{24.l3r+{^xlr% + 12/i^r^ + 15^0?'+ + 2) 

-13/x^r^ + ^0(24^14 + 22rt) + O^gr^ (IT^^r^ + 22) + I2filrl - 8) 

+l,lnrl{iilrl{m^r+ + 9f,lrl - 19) + 8^g(9/3r+ + 3^^r2 - 11) - 6^[]4 - 6fil) 

+2n'^r+{filrl{48^r+ + 3^^r2 - 25) + 3^^r2 (60/3r+ + Ufi^rl + 5) 

+/i^(72/3r+ + 27nlrl - 52) - 2//^r!|_ - 12/i^) + 320Vi^(12/3r+ + 7)) 

+0((r-r+)2), (C19) 



2 



J7r+ (0(4^2^2 _ g^2^2 ^ g) ^ r+(4^g + Sflj) + 403^2 + IQQ 



2 



6(r-r+)(^2^+ + 2f])2 
+ 6(»., + 1)(,2., I 2m,lr, + 40) (6^V(8/^r,(2.g.i + 3) - 5,2.2 ^ 4^2,2 _ 
+3//gJ7r2 (4/3//gr+ - fil) - 41)^2 (-24/3r+ + filrl + 18) 
+2!^3(i2/3r+(,^4 + 6,gr2 + 2) + 3/igr2 (,2^2 _ ^ g^2^2 _ ^g^ 

+02r-+(l2,^r2 (3;3r+ - 5) + ,o(48/3r+ + 3,ir2 - 82) - 12,^) 

+6n^rl-16n^rl) +0{r -r+). (C20) 
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